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Abstract: We investigate the dynamical effect of the cosmological constant A on a sin- 
gle spherical vacuum void evolving in the universe within a global solution of Newton- 
Friedmann models. As a result, the main characteristic is that the void expands with a 
huge initial burst up to match asymptotically the Hubble flow. The size of voids increases 
with f2 Q and with A, which interprets respectively by the gravitational attraction of borders 
from outside regions and by the gravitational repulsion of vacuum from the inner region. 
The A-effect on the kinematics intervenes significantly by amplifying the expansion rate at 
redshift z ~ 1.7 for a background density parameter Q ~ 0.3. For a class of parameters 
values, which corresponds in GR to spatially closed Friedmann models, it is interesting 
to note that a test particle in the inner region moves toward the border. Such a peculiar 
feature shows that the empty regions are swept out ; which stands as a stability criterion. 
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1. Introduction 

Foam like patterns with large empty regions are present in the distribution of observed 
galaxies at scales up to 100 Mpc (see e.g., |33], 16, |l^, [7], |26|, [| [l^, |j). Hence, one 



wonders whether the voids would be the most suitable structures for improving the study 
the formation process of large cosmological structures, and if their dynamics could be 
useful for determining the model of the universe. Investigations have been performed - 



on their statistical properties, by improving identification techniques (see e.g., [24]), by 
exploring their formation process in a ACDM model through N-body simulations (see e.g., 
fl3| , [32^ , 12, [2(], ||, p3|), and by probing their origins (see e.g., p2|| ); - on the kinematics 



of giant voids (see e.g., [|i8| ) ; - and on their dynamics, by testing models of void formation 
(see e.g., [|H|). 

At larger scales, according to observations, the dynamic of the cosmological expansion 
can be described by the Friedmann-Lemaitre (FL) model, despite the presence of large 
scale structures. The reason can be found on theoretical grounds, due to the fact that the 
recession of galaxies together with the black body spectrum of Cosmological Background 
Radiation ensure that FL model provides us with a suitable description of the universe 



Furthermore, the stability of FL models with respect to linear perturbations [19| enable us 



to use perturbative approaches in order to understand the large scale structures formation 
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process. Such investigations are generally performed with numerical simulations based on 
a FL cosmological background in which the structures are embedded. However, since voids 
are not linear structures [|J , a better understanding of the underlying dynamics requires a 
more adapted approach. With this in mind and the aim of understanding the effect of a 
cosmological constant A on the evolution of voids, we investigate the dynamics of a single 
empty spherical shell embeded in an uniformly distribution of dust which evolves with the 
cosmological expansion. Our model is based on solutions of Euler-Poisson equations system 
for whom kinematics satisfy Hubble (cosmological) law. Except for the description of the 
gravitational waves, FL models can be described within such a framework (see e.g. [pi]] ) and 
Newtonian treatment of perturbations leads to similar results as those obtained in general 
relativity (GR) Q. Analogous problems have been investigated in GR 
but with a different aim with respect to the cosmological constant effect on the dynamics 
of the void (empty bubble). 



27, 28 



3tt 31 



2. Euler-Poisson equations system 

The motion of a pressureless medium (dust) is described by its specific density p = p(f, t) 
and velocity v = v(f, t) at position r and time t. These fields are constrained by Euler- 
Poisson- Newton (EPN) equations system: i.e. Euler equations system 

^ + div(H = (2.1) 

dv dv _ , . 

where g = g(f, t) stands for the gravitational field, and they satisfy the modified Poisson- 
Newton equations system 

rot (f=0 (2.3) 
div<? = -AirGp + A (2.4) 

where G is Newton constant of gravitation. The sign of the cosmological constant is chosen 
(A > 0) such that the vacuum gravitational effect is repulsive. Accordingly to the status 
of gravitational field (which is not a vector field), we have chosen a cartesian coordinates 
system where g = at its origin, see |A[ 

In order to disentangle the evolution of structures from the cosmological expansion, 
it is convenient to write these equations with new coordinates, herein named reference 
coordinates. They are defined as follows 

(t, x = -), a = a(t) > 0, a(t) > 0, a Q = a(t ) = 1 (2.5) 
a 

where a is a solution of Friedmann type differential equation 

H + H*-± + ^p o a- 3 = 0, H = ± (2.6) 
3 3 a 
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and p stands for an arbitrary constant; herein, the dotted variables stand for time deriva- 
tives. This differential equation admits a integration constant 



K 



8ttG 



A 



3 - Po + --if o: 



H = H(a ) 



The integration of Eq. ( |2.6D provides us with Friedmann equation 



H 2 = - - ^ 



3 a 3 



> 0, 



lim H 

a— too 



(2.7) 
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where Hoc is the asymptotical value of the function H (a) when A > 0. 

The differential equation given in Eq. ( |2.8[ ) integrates by a quadrature, and the function 
a(t) is obtained as a reciprocal mapping. It depends on three parameters chosen among 



A, p , H Q and K D , see Eq. (|2.7| ). At this stage, it is important to mention that, excepted 
A, they do not identify with the cosmological parameters of FL model but they intervene 
solely in the choice of coordinates system. EPN equations system Eq. (ETT-2.4), reads in 
term of reference coordinates as follows 



+ div( / o c ?; c ) = 



dp c 
dt 

dv c dv c . 

-jj7 + ~^v c + 2Hv c 

at ox 



9c 



divg c 



rot g c = 

■-^ripc-Po) 



(2.9) 

(2.10) 
(2.11) 
(2.12) 



where "div" and "rot" denote the differential operators with respect to variable x, 



p c = pa , 



dx 
~dt 



(2.13) 



act respectively as the density and the velocity fields of the medium in the reference frame, 
and 

_ g /4ttG A^ 

Or. = \- " 



\3a* P °-3 ]X 



(2.14) 



as the acceleration field. 



2.1 Newton- Friedmann and Vacuum models 

In the reference frame, a uniform distribution of dust and an empty space are described, 
as obvious solutions of Eq. ( 2.9 , 2.10| , 2.11 . 2.12 ), respectively by : 



1. a Newton- Friedmann model (NF), which is defined by 

Pc = Po, v c = 6, g c = 



(2.15) 



Herein, according to Eq. (fT!. 2.13 ), a and H recover their usual interpretations as 
expansion parameter and Hubble parameter respectively, p Q = po? identifies to the 
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density of sources in a comoving space (of scalar curvature K a ). We focus our investi- 
gation to cosmological expansions a > 0, which excludes bouncing solutions. Hence, 
an analysis on roots of third degree polynomials shows that 

Kl < (4vrGp ) 2 A (2.16) 

must be fulfilled. The kinematics has two distinct behaviors, characterized by the 
sign of K Q . H decreases with time by reaching its asymptotic value either upward 
(if K Q < 0) or downward (if K a > 0) from a minimum defined by 



/ Kl 

H m = Hooi l - — — — - — < Hoc (2.17) 
V A(4vrG / 9 ) 

at a = AwGp K~ l , which defines the loitering period. Let us remind that FL models 
are characterized by the following dimensionless parameters 

Ao = ^, ^° = ^p k = ^ = X + n o -l (2.18) 

These notations are preferred to the usual £l\ = A and VLk = —k Q for avoiding 
ambiguities on the interpretation of cosmological parameters, see e.g., |8|| , []l0| . A 
dimensional analysis of Eq. (^7j) shows that the Newtonian interpretation of k cor- 
responds to a dimensionless binding energy for the universe during its expansion; the 
lower k the faster the cosmological expansion. In FL world model, k Q stands for 
the dimensionless spatial curvature of the comoving space; k D > characterizes a 
spatially closed universe. 

2. a Vacuum model (V), which is defined by 

Pc = 0, v c = {H^ ~ H) x, g c = ^f-p x ( 2 -19) 

3. Spherical voids in Newton- Friedman universe 

For avoiding misunderstandings in the interpretation of the dynamics, we use a covariant 
formulation of Euler-Poisson equations system (see ^]). The model of a spherical void 
embedded in a uniform distribution of dust is obtained by sticking together the V and NF 
models, both as local solutions of EPN equations system defined in Sec. |2.l| . It accounts for 
the dynamics of their common border (i.e. the boundaries conditions), which is assumed 
to be a shell of null thickness. 

3.1 Dynamical model 

We consider three distinct media : - a spherical thin shell (S) on which dust is uniformly 
distributed; - an empty space (V) inside; - and a uniform distribution of dust (NF) outside. 
For convenience in writing, each symbol S, V and NF denotes both the related medium 
and its model. They behave such that S makes the juncture of V with NF according to 
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Eq. ( [2.15|JO"9|) . The tension-stress on S is assumed to be negligible, which is characterised 
by a (symmetric contravariant) stress-energy tensor defined as follows 

T s 00 = (ps)c T^ = ( Ps ) c vi, Ti k = ( Ps ) c viv k c (3.1) 

The background is described by the following stress-energy tensor 

T™= Pc , T^ = 0, 2$ = (3.2) 

Because the eulerian functional 

T(x H+ 7 ) = J Tfj^dtdS + J TZl^dtdV (3.3) 

vanishes when 7 reads in the form 7^ = | (dp£ v + d u ^fA 1 one has 

+ gitj) + (djto + dotj - 2H& v{ + 4v k c d£ k ) ( Ps ) c dtx 2 dtt 

p c doto dtx 2 dxdn (3.4) 



'NF 

where dfi stands for the solid angle element (see 0). The radial symmetry of solutions 
enables us to write the reduced peculiar velocity and acceleration of a test particle located 
on the shell as follows 

v c = ax, g c = fix (3.5) 

where the functions a = a(t) and f3 = f3(t) must be determined. A by part integration of 
Eq. ( |3.4j) provides us with 

/ (<9o(ps)c + 3(ps)c« - Pcux) x 2 ^odt (3.6) 



t>2 



{do ((ps) c a) + 4(p s ) c a 2 + 2H(p s ) c a - (ps)^ 1 ) x 3 idt 



where x = \\x\\ stands for the radius of S and £ = + £2 + £2- This equality must be 
fulfilled for all bounded time interval and compact support 1-form. Hence, we easily derive 
the conservation equations for the mass 

do(Ps)c + (3(ps)c - p c x)a = (3.7) 

and for the momentum 

^+ (l + -^—x ) a 2 + 2Ha + - = (3.8) 



dt V G°s)c* 



With Eq. (2.5), the calculation of the gravitational force from the entire shell acting on a 
particular point 1 provides us with 



1 The modified newtonian gravitation field reads g = (~ 



A Gm\ - 
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It is noticeable that (ps)c = Pc x /3 ensures that the mass forming the shell comes from its 
interior 2 . It must be understood as a particular solution of Eq. ( p. 7] ) which is compatible 
with initial conditions provided by FL model. By focusing our investigation to such a 
model, Eq. (3.8) transforms 

da 9 27rG p r 

— +4a 2 + 2Ha ^ = (3.10 

at 3 a 6 

With the dimensionless variable 

X = 4^a 2 (3.11) 

where the ratio aH^ 1 stands for the expansion rate of S in the reference frame, Eq. ( 3.1C| ) 
tranforms into a Riccati equation 

(3 ' 12) 

where 

P(a) = A D a 4 - k Q a 2 + Sl a, P(l) = 1 (3.13) 
According to Eq. (p.l3| , |3.11| ), one has 

x = x,exp {£iSw (3I4) 

where X{ = x{ti) and = a(ij) stands for initial values (at time t = t{). 
3.2 Qualitative analysis 

A meaningful qualitative analysis of the expansion of the void must be performed by means 
of dimensionless quantities. Herein, we chose the magnification X and the expansion rate 
Y, which are defined as follows 

X = ~, Y = ^- (3.15) 

X{ ±1 

Their diagrams versus the expansion parameter a characterize adequately the expansion 
dynamics of S; they are obtained from Eq. ( [3.12 , 3.14) by numerical integration. They can 
be easily translated in term of cosmic time t since a(t) is a monotonic function of t in the 
present investigation. The initial conditions lie on the expansion rate Yi and the formation 



2 This property agrees with spherical voids originating from a uniformly distributed matter, although 
other behaviors can be envisaged. Indeed, Eq. (3.7) reads 

d S + 3a6 = 

where 

S = (ps)c - ^p c x = f(xi) exp ^-3 J adt 

accounts for the excess/defauft of the shell density with respect to that formed from an uniform distribution 
of matter, f(xt) stands for its initial value at t = ti. 



-6- 



date ti, as expressed by means of a, = a(tj). Herein, we investigate voids that initially 
expand with Hubble flow Yj = at a\ = 0.003. Instead of having an exhaustive analysis 
with respect to cosmological parameters, we limit our investigation around the generally 
accepted values A Q = 0.7 and Q = 0.3. 



The kinematics of the expansion is discussed in Sec. 3.2.1 , the dependence on cosmo- 
logical parameters (herein named Q-effect and A-effect respectively) in Sec. |3~2~2l , the role 



of initial conditions in Sec. 3.2.3, and the results are synthesised in Sec. 3.2.4. 



3.2.1 Kinematics 

The expansion velocity of the shell S with respect to its centre reads 

Y , H 
v = yHr, y = l + — h = — - 

h H 



(3.16) 



where y stands for the corrective factor to Hubble expansion. According to Fig. [l], which 
shows y versus a for Q = 0.3 and for three values of A € {0,0.7, 1.4}, S expands faster 
than Hubble expansion (i.e., y > 1) all along its evolution. It starts with a huge increase 
of the velocity at early stages (such as a burst) with no significant dependence on A, and 
then it decreases for reaching asymptotically the Hubble flow 3 . The A effect is to delay 




Figure 1: The corrective factor to Hubble expansion. 



the decreasing of the expansion velocity of S. On the diagram, it is characterised by a 
protuberance on the curves at z ~ 1.7 for A > 0.7 {i.e., k Q > 0) , the larger the A (i.e., 
k Q ) the higher the bump. Such a feature is due to the existence of a minimum value H m of 
Hubble parameter H, which is reached during the cosmological expansion and is related to 
a coasting period, see Eq. ( 2.17j ). With respect to observation of distant voids, the Doppler 



3 Although the present epoch (a = 1) appears quite peculiar because of the relative proximity of these 
curves, this turns out to be solely an artefact since they (will) cross at a > 1 but not in a single point. 
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shift between the edge and the center of a spherical void which is observed at redshift z 
does not exceed 

1 v c XY XjH Q 
Az=- -— = - — (3.17) 



[1 + z) c 1 + 2 



where c is the speed of the light. 



3.2.2 Dichotomy between cosmological parameters 

The growth of spherical voids is investigated by means of X and Y diagrams with respect 
to cosmological parameters o and A Q , with the aim of disentangling their related effects. 

1. O-effect. - The dependence on the outward density is analyzed with a vanish- 
ing cosmological constant A = and three values for the density parameter Q £ 
{0.15,0.3,0.45}. The comparison of diagrams in Fig.^ tell us that the larger the f2 Q 
the larger the magnification X. Such an effect results from the attraction of shell 



\ = o 




Figure 2: The magnification X — dependence on density. 

particles toward denser regions. According to Fig.||], the growth looks like a huge 
burst which freezes asymptotically up to matching Hubble expansion (Y = 0). More- 
over, the larger the Q Q the larger the expansion rate. This trend is not significantly 
modified by other acceptable values of A . 

2. A-effect. - The dependence on the cosmological constant is analyzed at constant 
S7 = 0.3 and A € {0, 0.7, 1.4}. The comparison of diagrams in Fig.|3| shows that the 
larger the A the larger the magnification X, which increases nonlinear ly with A . 
According to Fig. ||, the expansion rate Y does not characterize A since the related 
curves do not disentangle. Its dependence on A is weak and appears later on the 
evolution, but it has a cumulative effect which is reflected on X, see Fig. [I]. 

According to Eq. fl2.7|), the dependence of the dynamics of S on the parameter k Q can be 
retrieved from above analysis; the higher k the faster the S expansion, and a similar effect 
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Figure 3: The expansion rate Y — dependence on density. 
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Figure 4: The magnification X — dependence on the cosmological constant. 

for the magnification. Let us point out an interesting feature that characterizes the sign of 
k Q which appears on the motion in the reference frame of a test particle inside S. According 
to Eq. (EljjD : 

• if k Q < then the test particle moves toward the centre of the void; 

• if k a > then it starts moving toward its border at date t = t* defined such that 
a(t*) = ^nGpoK^ 1 , which happens when Hubble expansion H < H^. This feature, 
of sweeping out the inner part of S toward its border, stands for a stability criterion 
for voids which acts only for a spatially closed universe. 
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Figure 5: The expansion rate Y — dependence on the cosmological constant. 

3.2.3 Dependence on initial conditions 

- Figure H shows three curves corresponding to formation dates at € {0.003,0.03,0.3} by 
assuming Q a = 0.3, A = 0.7 and with Yj = 0. As expected, the earlier the date of birth 
the larger the magnification X. 



Q = 0.3 , X-7.7 




a 



Figure 6: The magnification X — dependence on the formation date. 

Figure^ suggests the existence of a limiting curve a \— > Y defined by ai = 0, which 
characterises £l - The expansion rate curves related to other formation dates are located in 
lower part and reach it asymptotically. One can expect that the effects on the dynamics of 
growth resulting from reasonable initial expansion rates 7^0, related to physical processes 
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Q = 0.3 , X-7.7 
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Figure 7: The expansion rate Y — dependence on the formation date. 



(e.g., supernovae explosions fig]), are negligible at primordial epoch 4 , what legitimises the 
initial condition Yi = at cij = 0.003. Evolutions with other initial conditions on Y^ can 
be deduced from that since any given point (aj,li) in the diagram belongs to a single 
evolution curve. 

3.2.4 Synthesis 

- As a result, the expansion of a spherical empty region does not show a linear regime, but 
it starts with a huge burst which freezes asymptotically up to match Hubble expansion. 
Its radius (in the reference frame) x increases from an initial size x% with the cosmological 
parameters f2 and A . The related individual effects can be interpreted respectively by 
the gravitational attraction from the outer parts and by the gravitational repulsion of 
its borders (or equivalently, the vacuum). The larger these parameters the higher the 
magnification X = x/x{. The dynamics is sensitive to f2-effect at primordial epochs and 
to A-effect later on. The evolution of its expansion velocity with time in the reference 
frame v c = YH Q x does not characterise A but £l . The cosmological constant intervenes 
mainly by delaying the decrease of the velocity expansion. The expansion rate evolutions 
related to formation dates show a common envelope curve which characterises the density 
parameter f2 D . 

3.3 Discussion 

It must be kept in mind that, based on two exact solutions of Euler-Poisson equations sys- 
tem, this model allows one to investigate properly the behaviour of a single spherical void. 
The connection conditions between the inner and the outer regions, which are investigated 
in a (covariant) Classical Mechanics, provides us unambiguously with the dynamics of the 



'Since Hubble expansion is all the more important towards the past, the earlier the formation date the 



weaker this effect, according to Eq. (3. IE) 
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shell. Such a schema does not correspond to an embedding of a void in a cosmological 
background solution (as usual), but it stands for a non linear approach of NF models. 
However, while it can be useful in N-bodies simulations to improve voids dynamics, it re- 
quires stringent properties on the space distribution of neighboring matter and ignores GR 
effects. At this step, one expects that it can apply solely to spherical voids with reasonable 
size in order to guarantee reliable results. 

With the aim of using the characteristics of voids to determine the values of the 
cosmological parameters, it is interesting to mention that the property of K a of being 
an integration constant (see Eq. 2.7 , 2.18| ) enables us to use k a to classify the dynamics 
of the cosmological expansion. With this in mind, let us note that the spatially closed 
world models k a > are characterized by a significant expansion of voids which reaches 
its maximum at redshift z ~ 1.7 (with S7 D ~ 0.3), the larger the A (i.e., k Q ) the higher 
the expansion rate. The inner part of voids shows a de Sitter expansion v = y 7 A/3r which 
sweeps out it, what interprets as a stability criterion. 



4. Conclusion 

Using Euler-Poisson equations system with cosmological constant, we investigate the dy- 
namics of an isotropic universe constituted by a spherical (empty) void surrounded by a 
uniform distribution of dust. It turns out that the behaviours of both regions (inside and 
outside the shell) coincide with Newton-Friedmann (NF) solutions, which ensures a consis- 
tent model. The dynamics of the void expansion is analyzed in the reference (comoving) 
frame 

The general behaviour of the void expansion shows a huge initial burst, which freezes 
asymptotically up to match Hubble expansion. While the corrective factor to Hubble law 
on the shell depends weakly on the cosmological constant A at early stages, it enables us to 
disentangle significantly cosmological models around redshift z ~ 1.7. The magnification 
of spherical voids increases with the cosmological parameters related to the density and 
to the cosmological constant. In case of spatially closed Friedmann models, the empty 
regions are swept out, which provides us with a stability criterion. It could be applied to 
observational data, as long as the matter distribution in the outer region of the single void 
is reasonably uniform 



A. Covariant formulation of Euler-Poisson equations system 



Souriau's covariant formulation of Euler-Poisson equations [|34| can be summarized as fol- 
lows : The geometrical interpretation of Newton dynamics || shows that the component 
of gravitational field g identify to the only non null Christoffel components of newtonian 
connexion 

1 oo _ 9 



Hence, one obtains their expression in the new coordinates system defined in Eq. ( |2.5D , and 
the only non null components read 

(r c )io = H5i (r c y 00 = 
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Let T be a function(al) defined on the set of symmetric density tensors x *— > 7 on newtonian 
spacetime R 4 

T(x^ 7 )= / T^-y^dtdV 

where didF stands for the volume element and T for a symmetric contravariant tensors 
which accounts for the media. The measure density T is eulerian if and only if it vanishes 
for all covariant tensor fields which reads 

lixv = j (duty + 

where d stands for the covariant derivative and x >— > £ for a compact support 1-form. In 
such a case, it is obvious to show that 

d„T^ = 

interprets as Euler equations, in any coordinates systems. 
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